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Abstract 

We study the classical and quantum motion of a relativistic charged 
particle on the spacetime produced by a global monopole. The self- 
potential , which is present in this spacetime, is considered as an ex- 
ternal electrostatic potential. We obtain classical orbits and quantum 
states for a spin- 1/2 and spin-0 particles. 
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1. Introduction 



Global monopoles are heavy objects probably formed in the early Universe 
by the phase transition that occcur in a system composed of a self-coupling 
scalar field triplet a whose its original global symmetry 0(3) is spontaneously 
broken to U(l). 

The simplest model that gives rise to a global monopole is described by 
the Lagrangian density below 

£ = l(d^ a ){d^ a ) - ^($ a $ a -iff (1) 

Coupling this matter field with the Einstein equation, Barriola and 
Vilenkin fl|] have shown that the effect produced by this object in the ge- 
ometry can be approximately represented by a solid angle deficit in (3 + 1)- 
dimensional spacetime, whose the line element of this manifold is given by 



ds 2 = -dt 2 + ^- + r 2 (d6 2 + sen 2 8d<p 2 ) , (2) 
a 

where the parameter a 2 = 1 — QitGrf is smaller than one and depends on the 
energy scale where the global symmetry is spontaneously broken. The area 
of a sphere of unit radius in this manifold is Aita 2 . 

Recently has been found that a charged particle placed in the spacetime 
produced by a global monopole becomes subjected to a repulsive electrostatic 
self-potential [0] given by 

V=« (3) 
r 

where r is the radial distance from the particle to the monopole and 

K= q ^±>0 (4) 

being S(a) a numerical factor which is positive for a < 1 and negative for 
a > 1. This self-potential is consequence of the topology of this spacetime 
and also its non- vanishing curvature. 

Here in this paper we decided to analyse the classical and quantum motion 
of a charged particle in the spacetime produced by a global monopole, taking 



2 



into account the self-potential as an external electrostatic potential in the 
equations of motion. So, with the objective to develop this analysis in a 
more general point of view we present this paper as follows: In section 2, 
we study the classical orbits by using the Hamilton- Jacobi formalism. We 
analyse the possibility of the test particle become bound to the monopole 
considering the unphysical case where a > 1. In section 3, we study the 
quantum motion considering the test particle as a bosonic and fermionic 
one, by using, respectively, the Klein-Gordon and Dirac equations. Finally 
in section 4 we summarize our main results and make some remarks. 



2. Classical Motion 



In this section we study the relativistic classical motion of a particle with 
mass M and electric charge q subjected to the electrostatic self-potential @ 
and in the manifold described by in (|2]). In order to do that, we shall 
use the Hamilton- Jacobi formalism. The first set of equations to be used is 
given by 

where 

A" = (6) 
qr 

being A a parameter along the trajectory of the particle. 
The second set is given by 

<r (H ~ ^) (If " ^) = ~ M2 - < 7 > 

For this system we shall employ the spherically symmetric ansatz for the 
functional S. 

S = -Et + R(r) + 6(0) + L z <p, (8) 

where R(r) and 0(8) are unknown functions and E and Lz constants of 
motion. 
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Substituting (||) into (0) we get the following set of differential equations: 

(de\ 2 li 



and 



9 fdR\ 2 ( „ iT\ 2 . L 2 , , 

a "* = - B + 7" M -^ (10) 



where L 2 is positive constant. 



We also can write 



d\ \ dr 
which, combining with (|l^) results in 



, r dr 9 / dR\ 



Also, from (|5]) we get 



M U * (13) 



and 



M I'D = < 14 > 



Combining ([13]) and (0) we obtain a differential equation relating the 
two angular variables, which by integration yields 



cot 2 6 = (x 2 ~ l)smV, (15) 
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with x — rj-We can see f rom the equation above that the relation between 
both angular variables does not depend on the parameter a. 

In order to analyze the dependence of the radial coordinate on the angular 
ones, what we call by equation of the classical orbits, we shall derive the 
function r((p). In order to simplify our calculation only, let us consider the 
surface 9 = tt/2, and defining u = l/r we get after some steps: 



a 2 x 2 




1 

L 2 



(E - Kuf 



AT 



whose solutions are 



u(<p) 



EK J(E 2 - M 2 )L 2 + M 2 K 2 



L 2 - 
for L 2 > K 2 , 



K 2 



u(<p) 



for L 2 < K 2 , and 



L 2 -K 2 



EK 



-cos 



y/L 2 - K 2 



a(<p - <fo) 



(16) 



(17) 



(E 2 - M 2 )L 2 + M 2 K 2 



K 2 -L 2 
cosh 



y/K 2 - L 2 

Tz 



K 2 -L 2 

a(<f - v^o) 



x 



u(<p) 



M 2 



Ea 2 x 2 (<p-cpoY 



(18) 



(19) 



for 



2EK 2K 

L 2 = K 2 . 

From the equations above we can see that the trajectories described by 
for the attractive case, K < ( which happens for the unphysical 
situation when a > 1), will be bounded from below by some positive value if 
E < M, so r will be bounded from above. If E > M, r can be infinity. The 
trajectories given by fllBD and (|T9D are spiral ones, and r goes to zero when 
if goes to infinity. In relativistic classical mechanics finite trajectories in 
general are not closed but rather rosetteshaped |J. The former case presents 
period given by 

2nL z 



T 



a 



VL 2 - K 2 ' 



(20) 
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which clearly on the parameter a. 

There is also a constant solution obtained by setting du/dip = in 
flTEp. This circular motion ocurs in relativistic theory just as in classical one. 

The relation between u = 1/r and the polar angle 9 can also be ob- 
tained. After some minor development from the previous equations we get 
for L 2 > K\ 



EK J(E 2 - M 2 )L 2 + M 2 K 2 

«M = -j]rr^ + jjtty, >< 



cos 



-a(ip - ipo) 



L z 

where the new angular variable is expressed in terms of 9 by 



(21) 




(22) 

We also should note that the equations above depend on the parameter 
a, however these two last relations differ from the previous ones envolving 
the azimuthal angle, where it appears as a factor multiplyng ip. 

Combining equation (|TTD with the fi = component of (|5|) we can obtain 
an expression for the radial velocity for the test particle. 

dr a 2 dR , . 

dt ~ E-K/rHr' ^ ' 

The turning points of the trajectory are given by setting dr/dt = 0. As 
consequence we get the potential curves given by 

—=d l2 y /2 — (24) 

M ~ \ + M 2 r 2 J + Mr' ^ ' 

For K > 0, which means the physical situation with a < 1 the right hand 
side of the above expression presents no extremals, so the particle cannot 
be bound by the monopole; however for K < 0, which is associeted with 
the unphysical situation with a > 1, it is possible to have E/M < 1 for a 
finite radius in (E5). (In this case the potential curves above coincides with 



6 



the circular trajectory mentioned previously). The situation where the self- 
potential is absent, K = , has been analysed in [Q. There it was found no 
possibility of the test particle be trapped by the monopole. 

After this analysis about the classical motion of the test particle in this 
manifold, where exact solutions were obtained, we shall leave for the next 
section the relativistic quantum motion study for this system. 

3. Quantum Motion 

This section is devoted to the relativistic quantum analysis of scattering 
and bound states of a massive charged test particle in the spacetime met- 
ric given by (E]), considering the self-potential as an external electrostatic 
potential. In this sense this section is a natural extension of the previous 
one, taking into account the quantum effects. In order to make this analysis 
complete, we shall consider the cases where the test particle is a bosonic and 
fermionic particle. We start first this analysis considering the former case. 

A) Bosonic Case 

The Klein-Gordon equation written in a covariant form in the presence 
of an external four-vector potential reads || 

[□ - l -^{d^g) - iq{d^) - 2iqA»d lt - q 2 A»A^ - M 2 ]y(x) = 0, (25) 

with 

Uipix) = -±=d,[V=-gg» v d v ip{x)l (26) 

being g = det(g^). 

We can also consider in our analysis the so called non-minimal coupling 
between the scalar field <f(x) with the geometry of the manifold itself. This 
term, which is invariant, is expressed by T R ip, where the coeficient T is 
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an arbitrary coupling constant and R the scalar curvature. (For the massless 
case the conformal coupling constant is T = 1/6.) 

The inclusion of the non-minimal coupling in our formalism, does not 
make our analysis more complicate. The reason for this fact is that the scalar 
curvature assaciated with the space-time produced by a global monopole is 
R = 2(1 — a 2 )/r 2 , which presents the similar behaviour of the centrifugal 
potential energy. 

Now, writing down the Klein-Gordon equation in the manifold described 
by (H), including the self-potential for a charged particle given by © and 
also the non- minimal coupling, we get from (|25|): 



a 2 L 2 K K 2 n 

[-d 2 + -,d r (r*d r ) - - - 2i-d t + -.-M 2 - T^Mx) = 0, (27) 



where 

7] = 2(1 -a 2 ) (28) 

and L the usual orbital angular momentum operator. 

Because our metric tensor is a static one and the self-potential is time 
independent, we shall adopt for the wave function the form below 

<p(x) = e- im R(r)Yr(e,<p), (29) 



where E is the energy of the particle. Substituting (p9|) into (|27|) we get 
the radial Klein-Gordon equation. 

H (^f ) - + ^ + ^-K 2 £- 2 -^R + (E 2 - M 2 )R = 0, (30) 



where we can identify the equation above with the Schrodinger one, associ- 
ating, 

T 2 L 2 + T V - K 2 
L -> { > 
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and 

E 2 — M 2 
~* 2Ma 2 ' 

We shall analyse in this section scattering states mainly; however if we 
admit the possibility of the parameter a 2 in (|2]) be bigger than one, the self- 
potential will be an attractive one, so bound states can be present in this 
case. Analysing the behaviour of the wave functions near the origin, and 
assuming E 2 > M 2 , well defined scattering states can be obtained expressed 
in terms of confluent hypergeometric functions as : 

R(r) = Cr^e^xFxiXi + l+i(3, 2(A, + 1); -2inr), (31) 



where A, = -\ + V^G+D + T,-g»] = = ek 



Following the standard procedure g we can take the asymptotic form of 
the confluent hypergeometric function obtain the long distance behaviour for 
the radial function. 



R ( ) r r(2A ' + 2)e ^ 2 C ° S{KT - (?r/2)(A; + 1} ~ f 31 ^ 2 ^ ± 111 (32) 
{T> ~ | T(A Z + l + i/3) \ nr ' 1 ' 



where ji = argT(\i + 1 + i(3). 

From the above equation it is possible to obtain the phase shift 5i, which 
is the most relevant parameter in the calculation of the scattering amplitude 



*i=|(l-A,)+7>- (33) 

As we can see from the result above, the phase shift presents two contri- 
butions : (i) From the modification of the effective angular quantum number 
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A; due to the geometry of the manifold itself and (ii) from the presence of 
the repulsive self-interaction term as another indirect consequence of this 
non-trivial topology. Taking K = rj = in ( |3"B| ) we reobtain the result found 
in the 0. 

Now, if we are inclined to consider the possibility of the test particle to 
be bound by the global monopole, we have to assume that the parameter 
a 2 in (||) is bigger than one. In this case the induced self-interaction is at- 
tractive and bound states can be obtained by taking in (^),i? 2 < M 2 and 
K = — | K | . Again imposing appropriate boundary condition on the solu- 
tions, they can be expressed as follows : 



R(r) = CX^-^F^Xi + 1 - £, 2(A ; + 1); 2nr) 



(34) 



where K = v^!EIZ and f = ££1. 

In order to obtain bound states we have to choose appropriately the 
parameters to terminate the series in fl3~4]). Admitting a polynomial of degree 
n for the hypergeometric function, we must impose the condition 



E\K\ 
a 2 K 



n. 



(35) 



With this condition we get discret values for the sel-energy given by : 



p(«) 



M 



a 2 (n + Xi + l) 2 
K 2 + a 2 (n + X l + l) 2 



1/2 



(36) 



with n = 0,1,2... 

Observing the expression above we call attention for two interesting re- 
sults: (z) The first one refers to the degeneracy problem. Because a^l, the 
self-energy depends on n + Az which is not in general an integer number, this 
reduce the degree of the degeneracy of our solutions, (ii) The second one 
refers with self- energy itself. As we can E^j depends on the parameter K 2 
which in general can be a very small quantity, so expanding (R3) in powers 
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of K 2 we get the non-relativistic limit of the energy [] 

K 2 



E n , ~ M 



2a 2 {n + \ l + l) 2 



(37) 



B) Fermionic Case 



The Dirac equation written in a covariant form and in presence of an 
external four-vector potential reads 



[i^{x){d a + iqA^ - - M]*(z) = 0, 



(38) 



where r M (x) is the spinorial connection coeficient given in terms of the tetrads 
and Christoffel symbols 



-)l [a h {h) eUd,e {b)v -T^e m ) 



(39) 



The generalized Dirac matrices 7 M (x) are given in terms of the standard 
flat spacetime gammas 7^ by the relation 



(40) 



For the metric corresponding to the global monopole we shall use the 
following tetrads : 



/ 1 















asinOcosip r 1 cos9cosip —(rsin9) 1 simp 



asinOsirup r cosOsinip 



\rsint 



COSif 



acosO 



r sinO 







(41) 



1 This result would be obtained if we have used the Schrodinger equation defined in the 
geometry defined in (Q), in the presence of an external electrostatic potential (|^). (See 
Ref. [2]) 



11 



which obey the relation 



e (a) e (b)V 



fj.li 



(42) 



being g fllJ (x) given by @. (Although there is a simpler basis tetrad to de- 
scribe the metric tensor associeted with a global monopole, the above choice 
is convenient, as we shall see, in the sense that taking a = 1 we reobtain the 
standard expression for the Dirac equation in a flat spacetime.) Now writting 
the Dirac equation in this spacetime and in presence of the self-potential we 
get: 



[i^) dt + i (r) dr+ W de + 

r rsmU 



1 <»>«.+ ^ 



T 



-q-f {0) A - M]^{x) = 0, 



(43) 



where 7^ = j.r,j^ = 7.6 1 and 7^' = j.cp, being r,9 and <p th standard 
unit vector along the three spatial directions in the spherical coordinates. 
We shall use the following representation of the flat space 7-matrices. 
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(0) 



1 
-1 



7 



(0 



a' 
-a 1 



and for the complete set solutions to the Dirac equation 

/ iF jm (r)<S> jm (6,<p) \ 



*(r,t) = - 
r 



iEt 



(44) 



which presents a well defined parity under the tranformation r 
The Qj m (8,(p) are the spinor spherical harmonics. 



— r. 



Substituting ([P|) into ( [13] ) and after some minor steps we find the set of 
radial diferential equations, 



(E — M — K/r)F, 



jm 



dGj m Gj r 

-a — / h r] 

dr r 



(45) 
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(E + M- K/r)G jm = + rj^f (46) 

where r\ = =)=(j + 1/2)0. Let us first consider the possibility of this system, 
charged particle and global monopole, presents bound states. This fact, as 
we have mentional previously, is only possible when K < 0, which implies 
a > 1. (Although this situations is realy an unphysical one, we decided to 
condider it below in order to make our analysis a complete one). We shall 
try solutions of the equations ( f45"l) and (|46|) in the form 

F jm (r) = C(l + E/M) 1 / 2 e- Kr (/ t r)^(F 1 (r) + F 2 (r)), (47) 

and 

G jm {r) = C(l - E/M)y 2 e- Kr (Kr) x ^{F 1 {r) - F 2 (r)), (48) 
where k = and Xj = VtHW-*' 

rv J rv 



Using in (|45|) and (f46"|) = — |X| and developing intermediate calcula- 
tions, we found that the solutions for the unknown functions Fi(r) and i^O") 
are expressed in terms of confluent hippergeometric functions as 



F 2 (r) = CiFi(A i -e,2A i + l;2Kr) 



and 



F 1 (r)=C 





-E\K\ 


M\K\ 


— rjdK 



1 F 1 (X j -£ + l,2X j + l;2 



KT 



(49) 



(50) 



where £ = EK/a 2 n. The expression for -F\(r) has been found using the 
general relations between the confluent hippergeometric functions II 

(z- — h a)iF 1 (a, b; z) = aiFi{a + 1, b; z). 
So the complete expressions to the functions F and G are given by : 



F jm (r) 



C(l + E/M)*e-" r (Kr) x * 



L Fi(Aj -e,2A i + l;2«r) + 



M|if| - T)OiK 



1^,--$ + l,2A,- + l;2/cr) 



(51) 
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and 



G jm (r) = -C(l - J E/M) 1 / 2 e- Kr ( K r)^j 1 F 1 (A J - £, 2X j + 1; 2«r) - 
a 2 «;A 1 ' — E\K\ , . I . 

- Wi-^ lFl(Aj ^ + 1 ' 2Aj + 1;2Kr) h (52) 

From the expressions above we can see the dependence of the radia func- 
tions on the parameter a which define this geometry and also on the external 
electrostatic potential through the constant K. 

In order to obtain discret values for the self-energies it is necessary to 
impose the vanishment condition on the wave function when nr^oo. From 
the asymptotic behaviour the confluent hippergeometric function M we get. 



E\K\ 



na 2 



\ j = n , n = 0,l,2... (53) 



Using the values for k and Xj given previously we obtain the explicit 
expression for the self-energy. 



E nij = M 



K 2 



(na + yftj + 1/2) 2 - K 2 f 



-1/2 

(54) 



Also from the expression above we see the energy depend on the parameter 
a in two different ways, through this parameter itself and the constant K 
indirectly. (Considering a = 1 and keeping the constant K ^ 0, the self- 
energy above reproduces the values of the self-energies of an electron in the 
hidrogen atom.) The explicit normalized eigenf unctions ([44]) can also be 
obtained. 

Let us now study scattering states. These states are obtained considering 
E 2 > M 2 in ( |45|) and (|46D . This is the real case because K > 0. 
Again we shall try solutions in the form : 

F j>m (r) = C(l + E/MY^iKr^iF^r) + F 2 (r)), (55) 
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G hm (r) = -C(l - E/M) 1 / 2 e iKr (Kr) x ^(F 1 (r) - F 2 (r)) 



(56) 



where « = «. 

a 

Substituting the expressions above in P5j) and PB]) and developing some 
intermediate calculations we find that i*i(r) and F 2 (r) can also be expressed 
in terms of complex confluent hippergeometric functions as : 



F 2 (r) = C 1 F 1 {X j + 2\j + 1; -2i«r) 



and 



Fi(r) = C 



a 2 (i/3 + Aj) 
+ z^a) 



-Mff 



i-Fi(Aj + l + i/3, 2Xj + 1; -2i«r), 



(57) 



(58) 



where /3 = EK/a 2 K. 

The complete expressions for the scattering wave function can be given 

by 



3,m 



C(l + E/M) 1/2 e iKr (/€r)^| 1 F 1 (A j + z/3, 2A^- + 1; -2«er) + 



a 2 (i/3 + Xj) 
+ irja) 



-MK 



xF^Xj + l + ip, 2Xj + 1; -2kr) , 



(59) 



and 



G 



G(l - ^/M) 1 / 2 e iKr (ftr) A ^|iF 1 (A i + ^,2Aj + 1; -2i«r) - 



a 2 (i/3 + Xj) 
irja) 



-MK 



iFx(Xj + l + 2Aj + 1; -2i«x) 



(60) 



As in the bosonic case we can also obtain the phase shift by the asymp- 
totic form of the wave function ([14]), which by its turn depend on the above 
expressions. So, let us write down the long distance behavior for F and G. 



F R-i 
- 1 jm ~ 



C (1 + E/M) 



i r(2A,- + 1) 



7T 



I COS 

[A, + 2/31 
2kt 



T(X J +t{3)\(-ii- l MK 

at - cr 2 



/37r/2 e |((7i-<T 2 ) x 



/«r — Aj— — [3ln(2Kr) + 7^ + ^ 



sm 



7T 



+ 

<J\ + cr 2 



Kr - A j'2 - /3/«(2/cr) + 7 3 - + o 



+ .4(61) 
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and 



G 



-C(l - E/M) 



i r(2A,- + 1) 



7T 



| sin 

|A 3 -+t/3| 
2«r 



J \ a m 2 



/37r/2 e i(<7i-<j 2 ) x 



nr — Xj— — f3ln(2Kr) + 7^ + 

_ _ 



cos 



Kr - Aj- - f3ln(2Kr) + jj + 



+ ... (62) 



where 



<T\ = arg 



(63) 



1 

° 2 = ar 9 7 ^HFi — 2 

and 

7j = ar(yfr(A J ' + (65) 
From the above expressions we can obtain the phase shift, which is given 

by 

^• = |[0' + l/2)-A i ]+7i + ^^- (66) 

As in bosonic case we can see that the phase shift presents contibuitions 
comming from : (i) The modification in the effective total angular quantum 
numbers, Xj, due to the geometry of the manifold it self and (ii) the presence 
of the induced self-interaction. There is also an additional contribution when 
we compare ( |63"D with (|33"D , the phase shift for the bosonic scattering states, 
this difference is due to specific behaviour of the spinor field in presence of 
an external potential. 



(64) 
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4. Concluding Remarks 



In this paper we have analysed the relativistic motion at classical and 
quantum point of view, of a charged particle in the neighbourhood of a 
global monopole, considering the induced electrostatic self- interaction as 
the zeroth component of an external four vector potential A^. Although the 
magnitude of the self-interaction can be a small quantity for a typical model 
of global monopole system comming from a grand unifield theory 0, which 
seems to be the more problabe, its effects can be measured by the phase shift 
for scattering elementary particle by a global monopole. 

In order to make our analysis as more complete as possible, we decided 
to investigate the behaviour of the classical and quantum motion of the 
charged particle, considering the possibility to have an attractive electro- 
static self-interaction. For this case we could observe that it is possible to 
construct bound states for the system composed by the particle and the 
global monopole. 

The main results of this paper were the obtainment of exact equation 
for classical orbits, self-energies, phase shifts and complete eigen-function at 
quantum level, for the movement of a charged particle in the neighbourhood 
of a global monopole. 

Finally we would like to make a comment about the procedure adoptel 
by us in this paper. Although we have used the induced electrostatic self- 
interaction given in Eq.(|3|) to study the relativistic motion a charged par- 
ticle, this result was obtained in Ref.[2] assuming that the particle was at 
a fixed position outside the global monopole. The complete procedure to 
determine the self-interaction due to a moving particle in the presence of a 
global monopole, should take into account the previous knowledgement of the 
particle's trajectory. Because this procedure makes a concret analysis about 
this system almost impossible we decided to adopt the quasistatic approach 
presented here, i.e, we neglected the corrections on the self-interaction com- 
ming from this movement Q 

2 For a typical grand unifield theory the parameter 77, associated with the scale where 
the global symmetry is spontaneously broken, is of the order 10 16 Gev and A = 1 — a 2 = 
8ttGt] 2 ~ 10~ 5 . In Ref. [2] it is estimated the magnitude of the induced self-interaction 

2 A 

for small A , it is approximately K w q * 2 . 

3 In Ref. [10] is presented the induced electromagnetic semmmlf-interaction on a charged 
particle moving along an arbitrary trajectory on a conical space-time. 
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